Coherent oscillations of current due to nuclear spins 
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We propose a mechanism for very slow coherent osciUations of current and nuclear spins in a 
quantum dot system, that may qualitatively explain some recent experimental observations. We 
concentrate on an experimentally relevant double dot setup where hyperfine interaction lifts the spin 
blockade. We study the dependence of the magnitude and period of the oscillations on magnetic 
field and anisotropy. 
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There are significant experimental and theoretical ef- 
forts aimed to utilize and manipulate the electron spin 
in the context of electronic transport, those are com- 
monly referred to as spintronics0,l3. Some of them in- 
volve nuclear spins as well. Many efforts concern GaAs 
semiconductor structures where the hyperfine interaction 
between electron and nuclear spin is relatively strong 
0. Furthermore, the nuclear spin relaxation times are 
much longer than the time scales related to electron 
dynamics^, 1^. This time scale separation has facili- 
tated experiments where a quasi-stationary polarization 
of the nuclear system was achieved and its effect on the 
electronic transport was observed 0, 0, 0|- Recent ex- 
periments implement novel ways of controlling the inter- 
action of electron and nuclear spin 0| , whereby the 
coherent oscillations between 'up' and 'down' polarized 
nuclear system have been observed . 

Transport experiments with quantum dots allow for a 
detailed study and control of the quantum dot energy 
spectrum, both in the regime of linear transport and 
in non-linear regime of excitation spectroscopy ^2 ll-^ - 
They reveal interesting regimes where unusual mecha- 
nisms of electron transport are the dominant ones. For 
instance, in the Coulomb blockade regime the direct 
charging of the quantum dot is forbidden by energy con- 
servation and the residual current is due to cotunneling 
[l^ ITsf . Another regime is the so-called spin blockade 
where the direct electron transfers are blocked by virtue 
of spin conservation 0| . Spin blockade may be achieved 
in various ways, e. g. with spin polarized leads p^. Re- 
cent experiment realizes the spin blockade in a double dot 
system, where the absence of transitions between spin- 
singlet and spin-triplet states in the dots explains the 
current rectification observed |l8|. Any spin-flip mech- 
anism facilitates these transitions, giving rise to a small 
residual current. 

The same group has recently reported [iflj l an unex- 
pected and unusual result. They have observed coherent 
oscillations of the residual current in this regime with a 
period in the range of seconds. This extremely long time 
scale together with the fact that the oscillation period 
and amplitude can be modified by resonant excitation of 
the nuclear spins, strongly suggests that the origin of the 



oscillations may be traced to the hyperfine interaction 

In this letter we propose a concrete mechanism for 
these oscillations that can at least qualitatively explain 
the observation made in Ref. The effect comes about 
from the dynamics of nuclear spins driven by hyper- 
fine interaction with electron spin. The configuration 
of the nuclear spins can be presented with two effective 
magnetic fields acting on the electron spin in the two 
dots. The difference of these fields lifts the spin blockade 
thereby affecting the average current and electron spin. 
These fields are quasistationary at the scale of succes- 
sive electron transfers. They precess around the external 
magnetic field (z-axis) with frequency 10^ Hz. Albeit this 
precession does not manifest itself in current oscillations. 
The oscillations result from slow nutations of the fields. 
These nutations arise from small deflection of electron 
spin from the z-axis, the deflection being induced by the 
fields. 

The dynamics of the nuclear spin fields appears to be 
far from chaotic so that the period and magnitude of 
the oscillations strongly depend on initial conditions of 
the nuclear spin system. Therefore the comparison with 
experiment may be only qualitative. In reality, the re- 
laxation of nuclear spin would lead to stabilization of the 
oscillations with a certain amplitude and period. How- 
ever, such stabilizing mechanisms would manifest itself 
at time scale much longer than the period. This is why 
we do not consider them in the present model. 

For details of the setup we refer the reader to Ref. 
In the regime of interest, the double dot can be 
in three distinct charge configurations as shown in Fig. 
n A charge configuration is characterized by {Nl^Nh), 
Nl{Nb) being the number of electrons in the left (right) 
dot. Transitions from (0,1) to (1,1) and from (0,2) to 
(0,1) are relatively fast involving electron tunneling ei- 
ther from the left or to the right lead with the rates Tlh. 
The bottleneck of the transport cycle are transitions be- 
tween (1,1) and (0,2). In the (0,2) configuration both 
electrons share the same orbital state, this makes it a 
non-degenerate spin singlet. In contrast to this, the (1,1) 
configuration comprises 4 possible states, grouped into 
spin-singlet and spin-triplet. The transitions between the 
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(1,1) singlet and (0,2) occur with the rate that is deter- 
mined by the tunnehng ampHtude between the dots and a 
relevant mechanism of inelastic scattering, e.g. phonons. 
These transitions do not require any spin- flip. The tran- 
sitions between the (1,1) triplet states and (0,2) singlet 
are forbidden by spin conservation: this causes the spin 
blockade. 

The hyperfine interaction with nuclei induces mixing of 
the singlet and triplet states in the (1,1) configuration. 
The part of the total Hamiltonian which contains the 
electron and nuclear spin operators reads 

H, = Kl-Sl + Kr ■ Sr + Ez ■ {SI + S'r) + 

Ast{Sl + SRf/2. (1) 

Here Sl,r are operators of electron spin in each dot [isj . 
The effect of nuclear spins is combined into two effective 
fields Kl R. In each dot 



Kl,R = f^'^^lL,R\kik] lL,R;f. 
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where Ik being operators of nuclear spin at Rk , the sum- 
mation goes over all nuclei and precession frequencies 'jk 
are set by an envelop of electron wave function and hy- 
perfine constant A. The third term in Eq. ^ represents 
the Zeeman splitting Ez in the external magnetic field 
while the fourth term represents the exchange splitting 
between the singlet and triplet. We adopt a semiclassical 
approximation of the effective fields Kl,r replacing them 
by classical variables^mii^. This approximation is 
justified by a big number of nuclei in the dots, iVqo ^ 1- 
The third and fourth terms in Eq. include the full spin 
only and therefore split the states onto singlet and three 
Zeeman-split triplet components. The mixing of these 
states is proportional to the difference of two effective 
fields Ka = Kl-Kr. 

The mixing thus lifts spin blockade. We assume that 
this mixing is the only mechanism of the residual cur- 
rent. This assumption is not crucial since alternative 
mechanisms, those include co-tunneling and non-nuclear 
spin-flips, would only produce an extra d.c. current back- 
ground for nuclear-induced current oscillations. 

We solve the problem in two steps. First, we solve for 
the density matrix of electron states assuming station- 
ary Kl,r- The output of the calculation are the average 
current and the average electron spin {Sl,r) in terms 
of Kl^r. Second, we use this output to derive equa- 
tions for the dynamics of K^ r and subsequently ana- 
lyze this dynamics. This approach relies on the time 
scale separation: the fields Kl r should not change at 
the time scale of successive electron transfers. The trans- 
fer rate can be estimated as {Ez — Ast) {Ka/ AstYTi. 
The small factor {Ka/Ast) is the ratio of the mixing 
and energy difference between singlet and triplet states 
and quantifies suppression of the current in the spin 



blockade regime. The fastest nuclear spin motion that 
changes K is the precession around external magnetic 
field with frequency lonmr — 10^ Hz. This results in 
the condition F^ {Ast/ KaY^nmr for the vahdity 
of our approach. As we see below, the current oscilla- 
tions are much slower with a typical period of the order 
of ')~^{A.sT / Ka)^ where 7 — {jl +7fl)/2. The preces- 
sion frequency and typical magnitude of effective field 
can be estimated ^ as 7 ~ En/NgD, K ~ En/ \/Nqd, 
where En « 0.135 meV in GaAs and A^qd ^lO*" for 
the quantum dots in question. The exchange splitting 
AsT oc lO^'"^ - 10~^eV as estimated in Ref. 19. This 
gives j~'^{Ast/Ka) ^ 0.01 - 0.1 sec. 

To make the first step, we describe the evolution of the 
electron system with Bloch equations for the density ma- 
trix approach. There are seven quantum states involved 
in the transport (|s) for the singlet in the (0,2) configura- 
tion, 1+) and |— ) for the two doublet components in the 
(1,0) configuration ,|0) for the singlet and |1), |2), |3) for 
triplet states in the (1,1) configuration, those correspond 
to = 1, 0, — 1 respectively) so that the full density ma- 
trix comprises of 49 elements. However, we can disregard 
most of the non-diagonal elements of the matrix except 
those between 4 states of the (1,1) configuration. So we 
end up with 19 equations only. We present here 3, this 
suffices to illustrate the overall structure. 
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where = {K% ± iK^)/V2. Note that the inelas- 
tic rate F; does not appear in Eq. (|3a|) . the same is 
true for the other triplets, but it is present in Eq. (|3b|) 
for the singlet. The average electron spins in each dot 
and the current can be readily obtained from the sta- 
tionary solution p"' of Eq. ©: {Sl,r) = Tt{p'''' Sl/r}, 
I = eF/fpl*. We do not need to present the cumbersome 
full solutions for the average spin here. Instead, we as- 
sume Tin ^ ^ L,R, AsT and note that in the zcroth 
order in Ka/Ast the average spin is in the z-direction 
and {xB = Ez/Ast) 
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With the same accuracy the current reads 
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In addition, we need the defiections of {Sl.r) from the 
z-direction. They arise in the next order in Ka/Ast and 



3 



are antiparallel in opposite dots: 




(6) 



Now we perform the second step and describe the dy- 
namics of the effective fields. To simphfy, we wih assume 
the same electron precession frequencies for all nuclei in 
each dot: jL,R-k = 1l,r- The advantage of this model is 
a closed system of equations for the collective fields: 



— Ka = Ja{Sa) X + 7GaAs-B X Ka 

at 



(7) 



{a — L, R) that describes precession of these fields, their 
moduli Kl u being constants of motion. The main pre- 
cession is around external magnetic field with the fre- 
quency 7GaAs|-B|- However, this precession is irrelevant 
since it does not change the if£ ^ and the angle (f> be- 
tween the projections of these vectors onto xy plane. 
These three slow variables actually determine the average 
spin and current, and the evolution equations for those 
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do not contain 7GaAs-B ^ 1l,r.{Slm) ■ These three equa- 
tions have an extra integral of motion, Kg = {Kf^'-^n + 
^rIl)/'^!- Moreover, the equation for two remaining 
variables appears to be of a Hamiltonian type. In di- 
mensionless variables k = K^/Kl,t = t^K^/ Ast the 
equations read 
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where 
£(fc,0) 
X{k) 
Y{k) 

m 



X{k) + Y{k) cos<^, 
l + kl- 2(fc|)2 



k + 2exBk^ + — , 



k^{l~{kl + k/2Y){kl-{kl-k/2)' 
{\ + 2xl)k^ + !^{l + kl-2{k^sf) 



Y{k) 
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and we introduced the notations k^^kg = 
Kr/Kl^K^Kl and e = (7^ - 7h)Ast/27^l. 
We also assume here that the asymmetry of electron 
precession frequencies is small, 5^ = |7l — 7i?,| ^ Il.r- 
One would expect this for the experiment in question 
since the two dots are nominally identical 

The C is evidently yet another constant of motion that 
depend on initial condition of the nuclear system. The 
solution of C = L, if it exists, determines a closed orbit in 



the (fc, (j>) phase space and the system experiences peri- 
odic motion along this orbit. This motion manifests itself 
in the current oscillations by virtue of Eq. |S1 The period 
and magnitude of the oscillations do depend on the initial 
conditions L,kii,kg. If L is close to 0, the period even 
diverges. There is, however, some regular dependence on 
the asymmetry e and the external magnetic field that en- 
ters through parameter x b , see Fig. |21 This dependence 
can be summarized as follows. The period T depends on 
asymmetry. If cxb ^1,7" — l""^ Ast I for a;^ < 1 and 
T ~ ^~^E\I AstH-l for xb ^ 1. In the opposite case of 
relatively large asymmetry tXB 3> 1, T ~ [xb^^)^^ for 
<C 1 and T ~ xb{^^)^^ for » 1. Also, the am- 
plitude of the oscillations relative to the average current 
is of the order of 1 for < 1 and decreases as x'^ for 
XB > 1. 

The initial values of K.a that determine the actual 
magnitude and period of the oscillations are distributed 
according to Gaussian statistics 'i^l with average squares 
corresponding to average squares of total nuclear spins in 
the dot. In Fig. |3|the current, see Eq. (Q, is plotted for 
various initial conditions but fixed e = 0.1 and xb = 1-6. 
We note apparent anharmonicity of the oscillations, this 
feature has been stressed in Ref. 

One might think that the periodic oscillations we ob- 
tained in our approach is an artefact of oversimplified 
model for nuclear dymanics in use. Recent work empha- 
sizes the importance of the fact that precession frequency 
7fe varies from nucleus to nucleus (2J]. This issue can be 
addressed within the semiclassical approach used here. 
To implement such approach numerically, we separate 
the nuclear spin system into iVf, ^ 1 blocks where the 
7fc 's are the same within each block but differ from block 
to block. The number of spins remains large so that the 
nuclear dynamics can be described in terms of effective 
fields ii-f^'^j^. This results in 6iVf, evolution equations sim- 
ilar to Eqs. IHl Intuitively, one expects such complicated 
dynamics to be chaotic, so that the memory about ini- 
tial conditions is lost after some time ~ 7^^. This would 
be really dreadful for the mechanism discussed, so we 
have performed extensive numerical simulations to check 
this circumstance [2^. To summarize the results, the 
dynamics is not chaotic, the memory about initial con- 
ditions persists and the nuclear system exhibits regular 
oscillations, typically with several periods. To illustrate 
this fact, we present a typical result in the inset of Fig. 
3. It shows the regular long-period motion and extra fast 
oscillations on the timescale of 7^^. 

In conclusion, we propose a mechanism whereby the 
transport via a double quantum dot induces slow regu- 
lar nutations of the nuclear spin system, these nutations 
are seen in the transport current. We model the con- 
crete experimental situation (|l9j|'l and our estimations 
of the typical frequency, anharmonic shape of the oscil- 
lation predicted, and sensitivity to magnetic field shown 
correspond to the observations made in [l9l |. More re- 
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search on relevant nuclear spin relaxation mechanisms 
is needed for detailed comparison with the experiment. 
From the other hand, the mechanism presented is suffi- 
ciently general and works in any conditions where the hy- 
perfine interaction provides the main mechanism of spin 
blockade lifting. 

We are grateful to the authors of Ref. |f9|J for commu- 
nicating their results prior to publication. We acknowl- 
edge the financial support by FOM. 
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FIG. 1: Each charge configuration is characterized by 
(NljNr) and also its multiplicity. The three charge trans- 
fer processes are (i) electron entering the left dot from the 
left lead, (ii) electron going from left dot to right dot, and 
(Hi) electron leaving the right dot to the right lead. See text 
for discussion on the F's. 
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FIG. 2: The period as a function the integral of motion L for 
various values of xb- The asymmetry is e — 0.1, ka = 0.9 
and kg — 0.0. The inset shows the period as a function of a;s 
for fixed values of L. 
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FIG. 3: The current shown as a function of time for e=0.1, 
xb=1-6, but for different initial conditions, i.e. different L. 
The inset shows current for L = —0.14 in the case of Ni,=l 
and 25, where Nb is the number of blocks. Note that the long 
period oscillations are still present for Nf, ^ 1. 



